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O ■ ABSTRACT 
(N ; 

c ■ 

^ I We consider a time dependent, spherically outflowing wind, in Schwarzschild 

^ I space-time, consisting of electron-positron pairs and photons . Without assuming 

^ I thermal equilibrium, we account for the microphysics, including two-body pro- 

^ I cesses (ee ee, ■ye ye, e~^e~ ^ 77) and their radiative three-body variants 

>■ I (ee ^ ee7, 7e ^ yey, e^e~ ^ 777)- We present a finite-difference scheme for 

^ ' solving the general relativistic kinetic Boltzmann equations for pairs and photons. 

^ ' We apply this to the concrete example of a wind from a hot, bare, strange star, 

■ predicted to be a powerful source of hard X-ray photons and e^ pairs created 

in : 

O ■ 

■ 

Qh' with total luminosities in the range 10'^^ — 10^^ ergs s ^, for different values of 

O ' the injected photon-to-pair ratio. The wind parameters-such as the mean op- 

' tical depth for photons, the rates of particle number and energy outflows, bulk 

■ velocity, and number density of the pair plasma-are presented as functions of the 
^ ■ distance from the stellar surface, as well as characteristics of the emergent radia- 

^ ' tion. We find that photons dominate in the emerging emission, and the emerging 

■ photon spectrum is rather hard and differs substantially from the thermal spec- 
trum expected from a neutron star with the same luminosity. This might help 
distinguish the putative bare strange stars from neutron stars. 

Subject headings: radiation mechanisms: thermal — plasmas — X-rays: stars — 
radiative transfer — stars: neutron 



by the Coulomb barrier at the quark surface. We study the kinetics of the wind 
particles and the emerging emission in photons and pairs for stationary winds 
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1. Introduction 

Compact astronomical objects identified with neutron stars, strange stars, and black 
holes are believed to be sources of electron-positron (e^) pairs that form in their vicinity by 
different mechanisms and flow away. Among these objects are radio pulsars (Sturrock 1971, 
Ruderman & Sutherland 1975, Arons 1981; Usov & Melrose 1996; Baring & Harding 2001), 
accretion-disc coronas of the Galactic X-ray binaries (White & Lightman 1989; Sunyaev et 
al. 1992), soft 7-ray repeaters (Thompson & Duncan 1995; Usov 2001b), cosmological 7-ray 
bursters (Eichler et al. 1989; Paczyhski 1990; Usov 1992), etc.. The estimated luminosity in 
pairs varies greatly depending on the object and the specific conditions: from ~ 10^^ — 10^^ 
ergs s~^ or less for radio pulsars up to ~ 10^° — 10^^ ergs s~^ in cosmological 7-ray bursters. 

For a wind out-flowing spherically from a surface of radius R there is a maximum 
(isotropic, unbeamed) pair luminosity beyond which the pairs annihilate signiflcantly before 
they escape (see Beloborodov 1999 and references therein). This is given by 

where F is the pair bulk Lorentz factor, and ctt the Thomson cross section. When the 
injected pair luminosity, L±, greatly exceeds this value the emerging pair luminosity , L± 
, cannot signiflcantly exceed L^^; in this case photons strongly dominate in the emerging 
emission: L± < L™^ <^ L± ~ L^. Injected pair luminosities typical of cosmological 7-ray 
bursts (e.g., Piran 2000), L± ~ 10^° — 10^^ ergs s~^, greatly exceed for their estimated 
F± ~ 10^. For such a powerful wind the pair density near the source is very high, and 
the out-flowing pairs and photons are nearly in thermal equilibrium almost up to the wind 
photosphere (e.g., Paczyhski 1990). The outflow process of such a wind may be described 
fairly well by relativistic hydrodynamics (Paczyhski 1986, 1990; Goodman 1986; Grimsrud 
& Wasserman 1998; Iwamoto & Takahara 2002). 

In contrast if L± <^ L™^^ annihilation of the outflowing pairs is negligible. It is now 
commonly accepted that the magnetospheres of radio pulsars contain such a very rarefled 
ultra-relativistic (F± ~ 10 — 10^) pair plasma that is practically coUisionless (see Melrose 
1995 for a review). 

Recently we described a numerical code and the results of calculations of spherically out- 
flowing, non-relativistic (F-j. ~ 1) pair winds with total luminosity in the range 10^^ — 10^^ 
ergs s"^ that is L ~ (10~^ - 10^)L™'''' (Aksenov, Milgrom, & Usov 2004, hereafter Paper I). 
(A brief account of the emerging emission from such a pair wind has been given by Aksenov 
Milgrom & Usov 2003.) While our numerical code can be more generally employed, the 
results we presented in Paper I were for a hot, bare, strange star as the a wind injection 
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source. Such stars are thought to be powerful sources of pairs created by the Coulomb barrier 
at the quark surface (Usov 1998, 2001a, see Fig. 1). For such luminosities it is not justified to 
assume thermal equilibrium, and so we have calculated the detailed microphysics including 
all two-body processes (ee ee, 76 — > 76, e+e" ^ 77) as well as their radiative three-body 
variants (ee ^ ee7, 7e ^ 7e7, e+e^ ^ 777)- The relativistic kinetic Boltzmann equations 
for pairs and photons were solved numerically to describe the time evolution and structure 
of pair winds and to find the emergent emission (we only presented results for stationary 
situations). Gravity of the star was neglected, and it was assumed that only pairs are 
emitted from the stellar surface. Thermal emission of photons from a bare quark surface, 
which has been neglected in Paper I, is strongly suppressed relative to black-body emission, 
but not completely, (e.g., Cheng &; Harko 2003; Jaikumar et al. 2004, sec Fig. 1). Still, 
whatever little is emitted may significantly affect the wind properties (see below). 

In the present paper we extend our previous studies to also include the effects of gravity 
by solving the Boltzmann equations for pairs and photons in a Schwarzschild space-time 
appropriate for the strange star. We also include thermal emission of photons from the 
stellar surface. 

In § 2 we formulate the equations that describe the pair wind and the boundary condi- 
tions. In § 3 we describe the computational method used to solve these equations. In § 4 we 
present the results of our numerical simulations. Finally, in § 5, we discuss our results and 
some potential astrophysical applications. 



2. Formulation of the problem 

We consider an pair wind that flows away from a hot, bare, unmagnetized, non- 
rotating, strange star. Space-time outside the star is described by Schwarzschild 's metric 
with the line element 

ds^ = -e'^f'c^df + e-^V + r\de^ + sin^ 9 dip^) , (2) 

where 

e'^ = (l-r,/r)V2 (3) 
is the lapse or redshift factor induced by gravity at a distance r from the stellar center, 

rg = 2GM/c^ ~ 2.95 x lO^M/M© cm (4) 

is the gravitational radius, with M the gravitational mass of the star. 
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FoUowing Page & Usov (2002) we consider, as a representative case, a strange star with 
M = IAMq which is constructed by solving the Tolman-Oppenheimer- Volkoff equation of 
hydrostatic equihbrium using an equation of state for SQM as described in Alcock, Farhi, 
& Ohnto (1986a) and in Haensel, Zdunik, & Schaefer (1986) [with a bag constant B = 
(140MeV)^, QCD couphng constant ac = 5'^/47r = 0.3, and the mass of strange quarks 
777.5 = 150 MeV]. The circumferential radius of the star is i? = 1.1 x 10^ cm. 

We assume that the temperature is constant over the surface of the star. The thermal 
emission in pairs and photons from the surface of strange quark matter (SQM) depends on 
the surface temperature, Ts, alone (Usov 1998, 2001a; Cheng & Harko 2003; Jaikumar et al. 
2004). The state of the plasma in the wind may be described by the distribution functions 
f±{p,r,t) and fj{p,r,t) for positrons (+), electrons (— ), and photons, respectively, where 
p is the momentum of particles. There is no emission of nuclei from the stellar surface, so 
the distribution functions of positrons and electrons are identical, /+ = /_. 

Since we fix the mass and radius of the central star the only remaining parameters that 
determine the wind are the surface temperature, which determines the pair injection rate 
and spectrum, and the photon surface radiation, which for default of better knowledge we 
take to be thermal with a suppression factor, rj, relative to black body; r) is our second 
free parameter. Because of the very steep dependence of the pair injection rate on the 
surface temperature, we actually use the former as a free parameter: in our calculations here 
L± = 10^^ 10^5, ... , 10^9 ergs s'K 



2.1. General Relativistic Boltzmann Equations 



The Boltzmann equations in a Schwarzschild background for either massless or massive 
particles can be obtained as a special case from the results of Harleston & Holcomb (1991) 
and Harleston & Vishniac (1992). These are given in a conservative form, which makes them 
particularly amenable to numerical treatment; 



d_ 
dfjL 



:i-/7v 



(5) 



Here, i is the species type {i — e for pairs and i = 7 for photons), ji — cos 9,9 \s the angle 
between the radius- vector r from the stellar center and the particle momentum p, p = |p|, 
13 f. — Vf./c = [1 — {nieC^ / eeYY^"^ ■, Pj — 1, and eg = c[p^ + imecf^-l'^ is the energy of electrons 
and positrons (for photons e-y = pc). Also, fd is the emission coefficient for the production 
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of a particle of type i via the physical process labelled by q, and Xi is the corresponding 
absorption coefficient. The summation runs over physical processes that involve a particle 
of type i. 



For convenience of numerical simulations we use instead of /j the quantities 

27re^f3ifi 



(6) 



standardly used in the "conservative" numerical method. This can provide exact conserva- 
tion of energy on a finite computational grid (see below). is the energy density in the 
{r, /X, e} phase space. Since all equations are the same for all particle species from now on 
we omit the index i. 

Prom equations (5) and (6) the Boltzmann equations can be written in terms of E as 



e--^ dE 


. d 


c at 


r^dr 








d/j, 







where 



E =J2iv'-x'E), 



(7) 



(8) 



2.2. Boundsiry conditions 

The computational boundaries are 

< r < Text, 0<t< t^t , (9) 

where Text = 1.66 x 10^ cm, and tst is the time when the wind approaches stationarity close 
enough, since here we concentrate on steady-state winds. 

Considering the interior of the quark star, it was shown that a thin (~ 10~^° cm) layer 
of electrons with an extremely strong electric field-the electrosphere-forms at the surface of 
strange quark matter (SQM) (Alcock et al. 1986a; Kettner et al. 1995), which prevents the 
electrons from escaping to infinity. The electric field in the elecrosphere is a few ten times 
higher than the critical field, Ecr ~ 1.3 x 10^^ V cm~^, at which the vacuum is unstable to 
creation of e~^e~ pairs (Schwinger 1951). Therefore, a hot, bare strange star is potentially 
a powerful source of e+e~ pairs created in the electrosphere and flowing away from the star 
(Usov 1998). We use in our simulations a pair injection rate (Usov 2001a): 

N'^ = 47ri?^F± , (10) 
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where 



F± = 3 X 10^^exp(-0.593C) 



X 



ln(l + 2C-^) , A 



+ 



;i + 0.074C)^ 2(13.9 + 



cm-V\ (11) 



and C = 20(Ts/10^K) ^. The energy spectrum of injected pairs is thermal with the surface 
temperature Tg, and their angular distribution is isotropic for < < 1. 

Thus, the pair injection luminosity, our first free parameter, is given by 

L± = 7V4-[mec2 + (3/2)A;Brs] , (12) 

where /cb is the Boltzmann constant (Usov 2001a). For the range we consider, L± — 10^^ — 
Iq39 gj-gg g-i^ ^jjg surface temperature is in a rather narrow range, Ts — (4 — 6) x 10^ K (see 
Fig. 1). 

Thermal emission of photons from the surface of a bare strange star is strongly sup- 
pressed if the surface temperature is not very high, Ts <S 10^^ K (Alcock et al. 1986a). The 
reason is that the plasma frequency of quarks in SQM is very large, fiWp ~ 20 — 25 MeV, 
and only hard photons with energies e-y > hwp can propagate in the SQM. The luminosity 
in such hard photons, which are in thermodynamic equilibrium with quarks, decreases very 
fast for Ts <^ huJ^/k-Q (Chmaj, Hacnscl, & Slominski 1991; Usov 2001a) and in our case, 
where Ts < 10^ K, it is negligible (sec Fig. 1). However, low-energy (e^ < hujp) photons may 
still leave the SQM if they are produced by nonequilibrium processes in the surface layer of 
thickness ~ c/cUp ~ 10~^^ cm (Chmaj et al. 1991). The emissivity of SQM in nonequihbrium 
quark-quark bremsstrahlung radiation has been estimated by Cheng & Harko (2003) who 
find that it is suppressed at least by a factor of 10^ in comparison with black-body emission, 
L-y ^ 10~^ Tbb- Usov (2004) and Usov, Harko, & Cheng (2005) have recently considered a 
modified model of the electrosphere taking into account surface effects for SQM . The mod- 
ification of the electrosphere results in an increase in the density of electrons by a factor of 
~ 20 — 30 in comparison with the case when the surface effects are ignored. The electron den- 
sity increase can additionally suppress the outgoing radiation from SQM in nonequilibrium 
quark-quark bremsstrahlung photons. In our simulations we take 

L^ = ?7Lbb = ?747rai?'TsS (13) 

where is a dimensionless free parameter. We present results with = 0, 3 x 10~®, 10~®. 
In default of better knowledge the photon spectrum is taken as black-body with the surface 
temperature Ts corresponding to L±. 
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An important parameter characterizing the affect of photon emission from the surface 
on the out-flowing wind is the Eddington luminosity -^Edd P^i^ plasma, above which 

the radiation pressure force dominates over gravity, 

L|dd ^ 10='^(M/1.4Mo) ergs s-\ (14) 

Within the narrow range of surface temperatures studied here, Ts ^ (4 — 6) x 10^ K, we get 
that is - (1 - 10)L|dd for = 3 x 10-^ and - (10^ - 10^)1^^^ for rj = 10"^ (see Fig. 1). 

The stellar surface is assumed to be a perfect mirror for both pairs and photons. At 
the external boundary (r = Text) pairs and photons escape freely from the studied region, 
i.e., the inward (// < 0) fluxes of both pairs and photons vanish there. 



2.3. Physical processes in the pair plasma 

As the plasma moves outwards photons are produced by pair annihilation {e'^e" — > 77). 
Other two-body processes that occur in the outflowing plasma, and are included in the 
simulations, are M0ller (e+e"*" — > e+e+ and e~e~ e~e~) and Bhaba {e'^e~ —>■ e~^e~) 
scattering, Compton scattering (76 — > je), and photon- photon pair production (77 e+e"). 

Two-body processes do not change the total number of particles in the system, and thus 
cannot in themselves lead to thermal equilibrium. So, we also include radiative processes 
(bremsstrahlung, double Compton scattering, and three-photon annihilation with their in- 
verse processes), even though their cross-sections are at least ~ ~ 10^ times smaller 
than those of the two-body processes {a — /he = 1/137 is the fine structure constant). 



3. The computational method 

Our grid in the {r, ji, e} phase-space is defined as follows. The r domain {R < r < Text) 
is divided into j'^ax spherical shells whose boundaries are designated with half integer indices. 
The j shell (1 < j < jmax) is between and rj+1/2, with Avj = r^+i/s - r^-i/a (ri/2 = R 

and rj^,,+i/2 = Text)- 

The //-grid is made of /c^ax intervals A^k = Atfe+1/2 - //fe-1/2: ^ <k < /Cmax- 

The energy grids for photons and electrons are both made of o^max energy intervals 
Ae^j = £u)+i/2 — £a;-i/2' 1 < <^ < <^max, but the lowcst energy for photons is 0, while that for 
pairs is mgC^. 
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The quantities we compute are the energy densities averaged over phase-space cells 

Eu,,kAt) - ^ f EdednT^dr. (15) 

where AX = Ae^A/ifcA(r|)/3 and A(r|) = r|_^i/2 - ^^-1/2- 

Replacing the space, angle, and energy derivatives in the Boltzmann equations (7) by 
finite differences (e.g., Mezzacappa & Bruenn 1993) we have the following set of ordinary 
differential equations (ODEs) for E^^^^j specified on the computational grid: 

^dE^ _^ A{r^fike^^E^^k)j 
c dt A(r|)/3 



where 



_ . 1 for photons, 

^'^ Ml-KcVe^)T/' for electrons, ^ ^ 



Ca;-l/2 + ^£^+1/2 
2 



lr=r 



(18) 



= (19) 

r, = 'Jrlllp^. (20) 

<A^ = <AW, .AH'A'W (21) 

Eoj,k{r) ^ — / E{e,i^,r)dedn, (23) 

E^Al^) = . \ 3 / ^(^' ^)drde, (24) 

Ek,j{^) ^ . \ 3 / E{e,ii,r)drdii, (25) 



A{r^l^ke^'^E^,k)j = r|+i/2e''^^+V2(;,^£;^^;,), 

-^|-l/2e''^^-^/^(/^fe^a,,fe)r-=r,_i/„ (26) 
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— ^U}-l/2Pui-l/2{Ek,j)uj-l/2, 



(27) 



A [(1 - fi')E^,j],^ (1 - /^Li/2)(^^J™+v2 

-l/2)(-^i^,i)M=Mfc-l/2 ' 



-(1 - lA-l/2){Eu,,j 



(28) 



' E. 



u}+l,k,j 

{^ui+2,k,j—^i^+l,k,3)i<^L^+l/2—^u>+l) 



{Ek,j)oj+l/2 - < 



li>0 

Eu),k,j 



4 

U<o, 



'w,k,j—J^ui-l,k,j)('^u>+l/2 



){^ui+l/2—^u>) 



(29) 



{Euj,j)fj.=^i^^i/2 — Euj,k,j 
^l^k{E^^k,j — E^^k-i,j) 
A^fc_i + A/Ik 



(30) 



{HkEuj^k)r 



"3+1/2 



(1 - XuJ,k,j+l/2) 



X 



Ew,k,j + 



E. 



+Xa;,fe,j+l/2/^fe- 



(31) 



Xw,A;,i+l/2 



1 + 



+ 



X(jJ,k, 



(32) 



The dimensionless coefficient x is introduced to describe correctly both optically thin and 
optically thick computational cells by means of a compromise between the high order method 
and the monotonic transport scheme without the artificial viscosity (e.g., Richtmeyer & 
Morton 1967; Mezzacappa & Bruenn 1993; Aksenov 1998). It is worth noting that if the 
values of (-E'fcj)w+i/2 and (£'(^j)^=^jfc+i/2 calculated by equations (29) and (30), respectively, 
are negative they are taken to be zero. 



The right side terms of equation (16) are 



(33) 
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UE)l,, = -^[ xE'^ded/.r'dr. (34) 

For the physical processes included in our simulations expressions rj^ ^ ■ and (x-B)^ ^ are 
calculated in Paper I. 

In paper I, where gravity is neglected, particle number is explicitly conserved when 
only two-body processes are taken into account; particle production occurs only via three- 
body processes. In the general relativistic scheme developed in this paper we were not able 
to achieve exact number conservation. However, in all runs for a stationary solution the 
change of the particle flux due to computational effects is at most a few percents. Energy 
conservation is still exactly satisfied. 

There are several characteristic times in the system. Some are related to particle reaction 
times, some to the time to reach steady state. These times may greatly differ from each 
other, especially at high luminosities (> 10^^ ergs s~^) when the pair wind is optically thick. 
The set of equations (16) is then stiff: at least some eigenvalues of the Jacobi matrix differ 
significantly from each other, and the real parts of the eigenvalues are negative. In contrast to 
Mezzacappa & Bruenn (1993) we use Gear's method (Hall & Watt 1976) to solve the ODEs 
finite differences analog (16) of the Boltzmann equations. This high-order, implicit method 
was developed especially to find a solution of stiff sets of ODE. To solve a system of hnear 
algebraic equations at any time step of Gear's method we use the cyclic reduction method 
(Mezzacappa & Bruenn 1993). The number of operations per time step is oc (cc;max/i;max)^imax, 
which increases rapidly with the increase of c^max and /Cmax- Therefore, the numbers of energy 
and angle intervals have to be rather limited in our simulations. 

Here we use an (e, yU, r)-grid with cjmax = 13, k^a.^ = 8, and jmax = 100. The discrete 
energies (in keV) of the e-grid minus the rest mass of the particles are 0, 2, 27, 111, 255, 353, 
436, 491, 511, 530, 585, 669, 766, and oo. This gives a denser grid at low energies and near 
the threshold of pair production, e = mgC^. The /i-grid is uniform, with AyU^ = 2/A;jnax = 1/4. 
The shell thicknesses are geometrically spaced: Ari = 2 x 10^'^ cm, and Ar^ = 1.3Arj_i. 

To check the effects of grid coarseness we also performed test computations with /Cmax = 4 
and 6, and separately with Wmax = 10 and 11. We did not observe changes in the results. 



4. Numerical results 

In this section, we give the results for the properties of spherically symmetric winds 
consisting of pairs and photons. The pair injection luminosity L± and rj are the parameters 
in our simulations. The photon injection luminosity is determined by equation (13) where 



-li- 



the surface temperature Ts relates to the value of L± [see equations (10) - (12)] . We start from 
an empty wind, injecting both pairs at a rate 10^^ ergs s~^ and photons at a corresponding 
rate. After a steady state is reached we start a new run with this steady state as initial 
condition, increase the energy injection rate in pairs by a factor of 10, wait for steady state, 
and so on. 

We next present the results for the structure of the stationary winds and their emergent 
emission at the external boundary. 

Figure 2 shows the mean optical depth T-y(r) for photons, from r to Text- The wind as a 
whole is optically thick [t^{R) > 1] for L± > 10^"^ ergs s~^ irrespective of 77. For Z± = 10^^ 
ergs s~^ the radius of the wind photosphere Tph, determined by condition T(rph) = 1, is 
~ 4 X 10^ cm (see Fig. 2). The wind photosphere is always deep inside our chosen external 
boundary (rph <^ ^ext), justifying our neglect of the inward (/x < 0) fluxes at r = Text- 

Figures 3 and 4 show the pair number density (rze) and the bulk velocity of the pair 
plasma onflow {v°^^), respectively, as functions of the distance from the stellar surface. For 
high luminosities {L± > 10^^ ergs s^^) both these quantities hardly depend on rj, i.e., the 
pair wind structure is determined only by L±. For low luminosities {L± < 10^"'' ergs s~^). 
He and Vg"* depend significantly on tj. In particular, for L± = 10^^ ergs s^^ near the stellar 
surface (r — i? < 10^ cm) the pair density for 77 = is ~ 6 times higher than for 77 = 10"^, 
while the velocity of the pair plasma is ~ 6 times smaller. This is because for a surface 
temperature Ts ~ 4 x 10^ K, which corresponds to L± — 10^^ ergs s~^, the mean kinetic 
energy of pairs near the surface, (3/2)A;bT's ~ 8.3 x 10~^ ergs is about half of the gravitational 
binding energy, so an atmosphere may form. In addition, for low surface photon luminosities 
(77 < a few X 10^^) radiation pressure is too weak to suppress the plasma atmosphere that 
results from the presence of gravity. So, an atmosphere does form, and, in turn, is conducive 
to pair annihilation. In contradistinction for rj = 10^^ radiation pressure does suppress the 
formation of an atmosphere, and the pair wind structure, i.e., the runs of rZe(r) and f °"*(r), 
is very similar to the case with no gravity and rj = 0. 

Figure 5 shows the rates of outflow of pairs {N±) and photons (N^/) through the 
surface as functions of r. For high luminosities {L± > 10'^^ ergs s^^) the rate of pair outflow 
N± almost doesn't depend on rj and decreases significantly at the distance 



from the stellar surface because of pair annihilation (see the estimate of Zann in Paper I) . For 




(35) 



low luminosities {L± < 10^^ ergs s ^) the dependence of N± on r is different for different 
values of 77. For 77 = the value of N± decreases outwards at least by a factor of 2 even when 
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L± is as small as 10^^ ergs s~^, while for rj — 10~^ the decrease of N± due to pair annihilation 
is very small for low luminosities, again reflecting the degree to which a pair atmosphere is 
formed. 

The photon outflow rate is seen to increase with increasing r not only because of pair 
annihilation, but also because of accumulative production of photons by radiative three-body 
processes. These processes are important for high luminosities (L-i- > 10^® ergs s"-*^) and are 
responsible for the increase of A^-y at r — i? ~ 10^ — 10^ cm for such high luminosities (see 
Fig. 5). 

The rates of energy outflow in pairs {E±) and photons (E^) vary with radius more 
or less similarly to the particle outflow rates, except that the total energy rate is explicitly 
conserved in all processes (see Fig. 6). The total energy rate decreases with increase of r 
because of the gravity effects, 

E(r) = E±{r) + E^{r) = e^^^-^)^ = \^^^L , (36) 

1 Tg/r 

and varies from L — L± + at r — R to the total emerging luminosity 

L = L± + = E{r^^) = ^^^^L (37) 

at r = Text- This differs qualitatively from Paper I where the gravity effects are ignored, and 
E{r) is constant. 

The number rates of emerging pairs (A^±) as functions of L± for different values of r] 
are shown in Figure 7. For L± = lO'^"^ ergs s~^ the value of N± for 77 = is ~ 3 times 
smaller than the same for rj = 10^^. The results of our new simulations for rj = 10^^ and 
our old simulations with rj = and no gravity practically coincide. This is consistent with 
our results on the wind structure. For L± > 10^^ ergs N± hardly depends on 77, and 
is ~ 1.5 — 2 times smaller than the result of Paper I (see Fig. 7). This is due to partial 
suppression of pair creation as the photon energies are reduced by gravitational redshift. 

Figure 8 shows the emerging total luminosities in pairs {L±, including the rest mass) 
and photons (L^)-given as fractions of the total emerging luminosity L — L± + L^-as func- 
tions of the injection luminosity, and for different values of 77. For the luminosity range we 
consider photons dominate in the emerging emission (L^ > L±), especially at high luminosi- 
ties where ^ L±. 

Figure 9 presents the energy spectra of the emerging photons for different values of 
L± and 77. For low luminosities {L± < 10^^ ergs s~^) and 77 = 0, photons, which form by 
pair annihilation, escape more or less freely from the star's vicinity, and the photon spectra 
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represent a very wide annihilation line that is redshifted in the gravitational field. The mean 

energy of the emerging photons varies from ~ 430 keV for L± = 10^^ ergs s~^ to 400 keV 
for L± = 3 X 10^^ ergs s~^ (see Fig. 10). If r] is near its maximum value (~ 10~^) there is 
no annihilation line in the spectra of emerging photons. Annihilation photons only modify 
these spectra by producing high-energy tails (see Fig. 9c). For high luminosities (L± > 10^^ 
ergs s^^) the energy spectra of emerging photons practically don't depend on r] and coincide 
with the results of Paper I. 



5. Discussion 

We have identified certain characteristics of pair winds outfiowing from hot, bare, strange 
stars and their emerging emission. For the energy injection rate in pairs we consider (L-t = 
10^"^ — 10'^^ ergs s^^) photons dominate in the emerging emission (see Fig. 8.). The spectrum 
of the emerging photons, we find, is rather hard (see Figs. 9 and 10) and differs qualitatively 
from the spectrum of the thermal emission from a neutron star with the same luminosity 
(e.g., Romani 1987; Shibanov et al. 1992; Rajagopal & Romani 1996; Page et al. 2004). In 
particular, for a bare strange star the mean energy of the emerging photons (see Fig. 10) 
is at least an order of magnitude larger than the same for a neutron star. This opens 
observational possibilities to distinguish strange stars from neutron stars. Hard spectra of 
bare strange stars are amenable to detection and study by sensitive, hard X-ray and soft 
7-ray instruments, such as INTEGRAL (e.g., Winkler et al. 2003). 

In this study we take into account both the thermal emission of photons from the strange 
star surface and gravity that have been neglected in Paper I. We have shown that for low 
values of ?7 < a few x 10~^ and L±. < 10'^^ ergs s^^ pairs emitted by the stellar surface are 
mainly captured by the gravitational field, and a pair atmosphere forms. The probability 
of pair annihilation increases because of the increase of the pair number density in the 
atmosphere, and this results in decrease of the fraction of pairs in the emerging emission in 
comparison with the case of Paper I (see Figs. 7 and 8). However, if the surface emission 
in photons is as high as its upper limit {rj ~ 10~^) radiation pressure forces dominate over 
gravity, and the pair wind structure practically coincides with the results for no gravity, 
irrespective of L±. 

The dominant physical processes in pair winds change significantly depending on L±. 
For L± < L^, 10^^ ergs the optical depth for photons is smaller than unity, r-y < 1 (see 
Fig. 2). In this case some part of e"^ pairs ejected from the stellar surface annihilates into 
photons in the process of their outflow, while photo-creation of pairs is rather rare. Above 
L± ~ pair creation by photons becomes important, and for L± ^ the rates of pair 
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creation and pair annihilation are more or less comparable at Zann ^ r — R <^ rph- At such 
distances from the surface the number density of pairs is nearly constant (see Fig. 3 and 
Fig. 8 in Paper I). In addition, for high luminosities {L^ > 10^^ ergs s~^) radiative three- 
body processes are important, and the total rate of the particle outflow increases with radius 
(see Fig. 5). These processes favour thermalization of pairs and photons in the outflow. 

For low luminosities (L± < L*) the spectrum of emerging photons significantly depends 
on the photo-emission from the surface, especially at low energies, <^ mgC^ (see Fig. 9). If 
surface photo-emission is negligible (rj = 0), this spectrum resembles a very wide annihilation 
line. The fractional emerging luminosity in the annihilation line decreases with the increase 
of rj. For rj = 10~^ the annihilation line is completely washed out, and a high-energy part 
of this line is observed as a high-energy tail (see Fig. 9c). For high luminosities {L± > 10^^ 
ergs s~^) the spectrum of emerging photons is practically independent of the photon emission 
from the stellar surface. 

Soft 7-ray repeaters (SGRs), which are the sources of short bursts of hard X-rays with 
super- Eddington luminosities (up to ~ 10^^ — 10^'^ ergs s~^), are potential candidates for 
strange stars (e.g., Alcock et al. 1986b; Cheng & Dai 1998; Usov 2001c; Ouyed et al. 2005). 
The bursting activity of SGRs may be explained by fast heating of the stellar surface up to 
a temperature of ~ (1 — 3) x 10^ K and its subsequent thermal emission (Usov 200 lb, c). 
The heating mechanism may be either impacts of comets onto bare strange stars (Zhang 
et al. 2000; Usov 2001b) or fast decay of superstrong (~ 10^"^ — 10^^ G) magnetic fields 
(Thompson & Duncan 1995; Heyl & Kulkarni 1998). Since for high luminosities (> 10^^ 
ergs s^^) the spectrum of emerging photons practically doesn't depend on the surface photo- 
emission, which is poorly known, it is possible to calculate more securely many properties 
(hght curves, energy spectra, etc.) of powerful (L » 10^^ ergs s~^) X-ray bursts expected 
from the stellar surface heating neglecting surface photo-emission. We plan to apply the 
tools developed here to study these problems. 

We thank members of the Israeli Center for High Energy Astrophysics, D. Eichler, A. 
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the Israel Science Foundation of the Israel Academy of Sciences and Humanities. 
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Fig. 1. — Injection luminosities of a hot, bare, strange star in e+e^ pairs (dotted line), 
in thermal equihbrium photons (dashed line), and the total (sohd line) as functions of the 
surface temperature Tg. The theoretical upper limit on the luminosity in non-equilibrium 
photons, 10~^Lbb (Cheng & Harko 2003), is shown by the dot-dashed line, Lbb being the 
blackbody luminosity. 
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Fig. 2. — Mean optical depth for photons, from r to infinity, as a function of the distance 
from the stellar surface for 77 = (solid lines), 77 = 3 x 10~^ (dashed lines), and r) — 10~^ 
(dotted lines), and for different values of the injected pair luminosity, L^, which increases 
in steps of factor ten from 10^^ ergs s~^ (lowest triplet of lines) to 10^^ ergs s~^ (uppermost 
triplet of practically coinciding lines) . 
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Fig. 3. — Pair number density as a function of the distance from the stellar surface, hne 
designation as in Fig. 2. 
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Fig. 4. — Bulk velocity of the pair plasma as a function of the distance from the stellar 
surface for (a) 1] = 0, (b) 77 = 3 x 10~^, and (c) rj = 10~^, shown for different values of L±, 
as marked on the lines. 
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Fig. 5. — Particle number outflow rates in photons (solid lines) and pairs (dashed hues) 
as functions of the distance from the stellar surface for (a) 77 = 0, (b) 77 = 3 x 10~^, and (c) 
rj — 10~^, shown for different values of L±, which increases in steps of factor ten from 10^^ 
ergs s~^ (lowest hues for each species) to 10^^ ergs s""^ (uppermost hues). 
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Fig. 6. — Rates of energy outflow, as in Fig. 5. 
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Fig. 7. — Number rates of emerging pairs as functions of the injected pair luminosity for 
77 = (thin sohd hne), 77 = 3 x 10~^ (dashed hne), 77 = 10"® (dotted hne). The result of 
Paper I where gravity has been neglected and 77 = is shown by the thick solid line. 
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Fig. 8. — Fractional emerging luminosities in photons (upper triplet of thin lines) and pairs 
(lower triplet of thin lines) as functions of the injected pair luminosity for = (thin solid 
lines), ?7 = 3 X 10~^ (dashed lines), and r] = 10'^ (dotted lines). The results for the case 
where gravity is neglected and 77 = are shown for comparison by thick solid line and thick 
dot-dashed line for emerging photons and emerging pairs, respectively. 
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Fig. 9. — Energy spectrum of emerging photons for (a) rj = 0, (h) r] = 3 x 10^^, and (c) 
1] = 10"^, shown for different values of L±, which increases in steps of factor ten from 10^^ 
ergs s~^ (lowest lines) to 10^^ ergs s~^ (uppermost lines). 
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Fig. 10. — Mean energy of the emerging photons (thin sohd hnes) and electrons (thin dashed 
hnes) as functions of the injected pair luminosity for different values of 77, as marked on the 
lines. The results of Paper I where 77 = and no gravity are shown by the thick solid and 
thick dashed lines for photons and electrons, respectively. 



